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D ■ Abstract 

We study the ground-state properties of the 5 = 1/2 antiferromagnetic Heisenberg model on 
the Union Jack strip lattice by using the exact-diagonalization and density matrix renormaliza- 
tion group methods. We confirm a region of the intermediate-magnetization state between the 
Neel-like spin liquid state and the conventional ferrimagnetic state of Lieb-Mattis type. In the 
intermediate-state, we find that the spontaneous magnetization changes gradually with respect to 
the strength of the inner interaction. In addition, the local magnetization clearly shows an in- 
commensurate modulation with long-distance periodicity in the intermediate-magnetization state. 
These characteristic behaviors lead to the conclusion that the intermediate-magnetization state is 
the non-Lieb-Mattis ferrimagnetic one. We also discuss the relationship between the ground-state 
properties of the S = 1/2 antiferromagnetic Heisenberg model on the original Union Jack lattice 
and those on our strip lattice. 
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I. INTRODUCTION 



Ferrimagnetism is a fundamental phenomenon in the field of magnetism. The most 
famous type of ferrimagnetism is called Lieb-Mattis (LM) onelH6[. For example, this fer- 
rimagnetism appears in the ground state of the (s, S)=(l/2, 1) mixed spin chain with 
nearest-neighbor antiferromagnetic interaction. In this system, the occurrence of the LM 
ferrimagnetism originates from the situation that two different spins are arranged alter- 
nately in a line owing to the AF interaction. In the LM ferrimagnetic state, the spontaneous 
magnetization occurs and the magnitude is fixed to a simple fraction of the saturated mag- 
netization. As in the case of this mixed spin chain, not only the magnetic properties but 
also the occurrence mechanism of the LM ferrimagnetism are well known since this type of 
ferrimagnetism has been studied extensively. Especially, the ferrimagnetism in the quantum 
Heisenberg spin model on the bipartite lattice without frustration is well understood within 



the Marshall-Lieb-Mattis (MLM) theorem 



On the other hand, a new type of ferrimagnetism that is clearly different from the LM 
ferrimagnetism has been found in the ground state of several one-dimensional frustrated 
Heisenberg spin systems [7- 13]. The spontaneous magnetization in this new type of ferri- 
magnetism changes gradually with respect to the strength of frustration. In addition, the 
incommensurate modulation with long-distance periodicity in local magnetizations is ob- 
served as a characteristic quantum behavior of the new type of ferrimagnetism. Hereafter, 
we call the new type of ferrimagnetism non-Lieb-Mattis (NLM) type. The mechanism of the 
occurrence of the NLM ferrimagnetism have not yet been clarified in contrast to the case of 
the LM ferrimagnetism. 

Historically, some candidates of the NLM ferrimagnetism among the 2D systems were 
already reported. For examples, there are the mixed-spin J1-J2 Heisenberg model on the 



square lattice 14[ and the S = 1/2 Heisenberg model on the Union Jack lattice of Fig. []Ja) 
15KL8fl. These 2D frustrated systems have the intermediate ground-state, namely "canted- 
ferrimagnetic state" as described in Fig. [2} in which the spontaneous magnetization is 
changed when the inner interaction of the system is varied. It has not been, however, in- 
vestigated whether the incommensurate modulation with long-distance periodicity exists or 
not in the local magnetization of the intermediate-magnetization state owing to the diffi- 
culty of treating these 2D frustrated systems numerically and theoretically. Therefore, the 
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FIG. 1: (Color online) Structures of the lattices: the Union Jack lattice (a), the Union Jack strip 
lattice (b). An S = 1/2 spin is located at each site denoted by a black circle. Antiferromagnetic 
bonds J\ (bold straight line) and J2 (dashed line) are represented. Sublattices in a unit cell of 
lattice (b) are represented by A, A', B, C, C, and D. 

relationships between the intermediate-magnetization states of these 2D frustrated systems 
and the NLM ferrimagnetic state are still unclear. 

Under such circumstances, quite recently, the S = 1/2 antiferromagnetic Heisenberg 
model on the spatially anisotropic kagome lattice was studied[19(. In this model, the 
intermediate- magnet ization states exist between the LM ferrimagnetic state and the non- 
magnetic one 



20 



24j . It was reported that the local magnetization in these intermediate- 



state shows large dependence on the position of the sites although it is difficult to judge 
clearly whether the incommensurate modulation with long-distance periodicity is present or 



absent. In addition, the S 
strip lattices were studied 
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/2 Heisenberg models on the quasi-one-dimensional kagome 



26|. These strip lattices share the same lattice structure in 



their inner part with the spatially anisotropic kagome lattice. The local magnetizations in 
the intermediate-state clearly show incommensurate modulations with long-distance peri- 
odicity irrespective of the strip width. Therefore, these results strongly suggest that the 
intermediate-magnetization states not only of the kagome strip lattices but also of the orig- 
inal kagome lattice are the NLM ferrimagnetism. 

These kagome results motivate us to investigate the ground-state properties of the quasi- 
one-dimensional strip model whose lattice structure is common to the part of the 2D lattice 
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FIG. 2: (Color online) Ground-state phase diagram of the 5 = 1/2 Heisenberg model on the 
Union Jack lattice depicted Fig. Ufa). Here, black and white circles represent up-spin and down- 
spin respectively. 



known as the other candidates of the NLM ferrimagnetism. In this study, we treat the 
5 = 1/2 Heisenberg model on the Union Jack strip lattice depicted in Fig. QJb). This strip 
lattice share the same lattice structure in the inner part with the original Union Jack lattice 
depicted in Fig. (Ha). Our numerical calculations lead to the conclusion that the NLM 
ferrimagnetic phase appears in the ground-state of the Union Jack strip model of Fig. [U(b). 
We also discuss the relationship between the ground-state properties of the present strip 
model and those of the original 2D model. 



II. MODEL 

The Hamiltonian of the 5 = 1/2 antiferromagnetic Heisenberg model on the Union Jack 
strip lattice depicted in Fig. [TJ^b) is given by 

% = J\ / JS^A ■ Sj 5 c + Sj_B • Sj 5 D + Sj t A' • Si,C' + Sj,A • S^b + S^b • Sj,A' 

i 

+ &i,C • Sj t D + S-^d • S^c + S^c • Sj+i^A + S^d • Sj+^B + S^c • Sj + i 5 A'] 

+ Jl ^^[Sj,A • Sj 5 D + S^A' • S^d + S^d • Sj+i^ + S^d • Sj+^A'], (1) 

where is an 5 = 1/2 spin operator at £-sublattice site in z-th unit cell. Positions 
of the six sublattices in a unit cell are denoted by A, A', B, C, C and D in Fig. QJb). 
We fixed J\ — 1 hereafter as a energy scale. In what follows, we examine the region of 
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< Jil J\ < 3.5 in the present study. Note that the number of total spin sites is denoted 
by iV; thus, the number of unit cells is N/6. 

Let us introduce here the ground-state phase diagram of the S = 1/2 Heisenberg model 
on the original Union Jack lattice of Fig. (U^a). At small J2/J1, one can see immediately that 
the antiferromagnetic Neel order is observed since this model corresponds to the S = 1/2 
antiferromagnetic Heisenberg one on the simple square lattice in the limit of J2/J1 = 0. 
When the J2/J1 is increased, the intermediate- magnetization state appears. A variational 
analysis for the classical model revealed that the spin configurations in this intermediate 
region are defined by the angle of cant (p as illustrated in Fig. I2"jl8 |. Therefore, this 



intermediate-state is called canted-ferrimagnetic one. The phase transition between the Neel 
and canted-ferrimagnetic phases occurs at ot\ = J2/J1 ~ 0.84 from the view point of the spin- 



wave theory 15|, 1 161 ] . On the other hand, it was reported that this phase transition occurs 
at a\ ~ 0.65 by using the series expansion (SE)0 and coupled cluster method (CCM)Q| 



techniques. It was also discussed the possibility that the semistriped-ferrimagnetic state as 



illustrated in Fig. |2] appears at very large values of J2/J1 («2 = J2/J1 ~ 125) 18] . 

In what follows, we examine the ground-state phase diagram of the 5 = 1/2 Heisenberg 
model on the Union Jack strip lattice depicted in Fig. H^b) and compare the results of the 
original 2D lattice with those of the our strip lattice. 

III. NUMERICAL METHODS 

We employ two reliable numerical methods: the exact diagonalization (ED) method and 



the density matrix renormalization group (DMRG) method 271. |28fl - The ED method can 
be used to obtain precise physical quantities for finite-size clusters. This method does not 
suffer from the limitation of the shape of the clusters. It is applicable even to systems with 
frustration, in contrast to the quantum Monte Carlo (QMC) method coming across the 
so-called negative-sign problem for systems with frustration. The disadvantage of the ED 
method is the limitation that available system sizes are very small. Thus, we should pay 
careful attention to finite-size effects in quantities obtained from this method. 

On the other hand, the DMRG method is very powerful when a system is (quasi-)one- 
dimensional under the open-boundary condition. The method can treat much larger systems 
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FIG. 3: (Color online) (a) Dependences of the lowest energy on S* ot . Results of J2/J1 = 0.1, 1.5 
and 3.5 for the system size of ./V = 48 are presented. Arrows indicate the values of the spontaneous 
magnetization M in each 3%j J\. (b) J2/ Ji-dependence of M/M s obtained from ED calculations for 
N = 24 (black cross) under the periodic-boundary condition and DMRG calculations for TV = 48 
(red triangle), 96 (blue square) and 144 (green pentagon) under the open-boundary condition. 

than the ED method. Note that the applicability of the DMRG method is irrespective of 
whether or not systems include frustrations. In the present research, we use the "finite- 
system" DMRG method. Note that we carefully choose the maximum number of retained 
states (MS) and the number of sweeps (SW) in our DMRG calculations. 

IV. RESULTS 

In this section, we present our numerical results in the ground-state of the S — 1/2 
Heisenberg model on the Union Jack strip lattice of Fig. \Bjo). First, let us explain the way 
to obtain the spontaneous magnetization M in the ground state of the quantum system with 
isotropic interactions. We calculate the lowest energy E(J 2 /Ji, S* ot , N), where S% ot is the 
z-component of the total spin. For example, the energies for each S* ot in the three cases of 
J2/J1 are shown in Fig. Ufa). In this figure, the results of the DMRG calculations with the 
MS = 700 and SW = 15 are presented when the systems size is N = 48 for J 2 /Ji=0.1, 1.5, 
3.5. The spontaneous magnetization M(J%/ J±, N) is determined as the highest «Sf ot among 
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FIG. 4: (Color online) (a) Spin configuration in the Neel-like spin liquid phase, (b) Spin configu- 
ration in the LM ferrimagnetic phase of M/M s = 1/3, where this configuration is obtained from 
the numerical results of the local magnetization shown in Fig. [5(b). 

those at the lowest common energy [see arrows in Fig. [3(a)]. 

Our results of the J2 / Ji-dependence of the M/M s are shown in Fig. [3(b), where M s 
means saturated magnetization value, namely, M s = N/2. In the limit of J2/J1 = 0, this 
Union Jack strip model depicted in Fig. QJb) is reduced to the 5 = 1/2 antiferromagnetic 
Heisenberg model on the three-leg ladder as known the typical system of the gappless spin- 



liquid ground states 
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301 ] . According to the study of the 5 = 1/2 frustrated three- leg 



spin ladder [31[, it is expected that the Neel-like spin liquid phase occurs in the ground-state 
when the strength of the J2/J1 is small but finite, where the schematic spin configuration in 
the Neel-like spin liquid state is illustrated in Fig. Ufa). Indeed, our numerical calculations 
lead to the same conclusions that the nonmagnetic phase of M/M s = appears in the region 
where J2/J1 is relatively small. 

For larger J2/J1, on the other hand, the magnetic phases with M/M s 7^ appears in 
the ground state. Careful observation enables us to find that there are two magnetic phases 
in the thermodynamic limit; one is the intermediate magnetic phase of < M/M s < 1/3 
and the other is the phase of M/M s = 1/3. It should be noted here that the phase of 
M/M s = (| — jj) which is found only under the open-boundary condition merges with the 
phase of M/M s = 1/3 in the thermodynamic limit of iV — > 00 since the value of M/M s 
becomes gradually larger and approaches the value of M/M s = 1/3. This change due to the 
increase of the system size comes from finite-size effect. It is important that we successfully 
observe the intermediate- magnetization phase where the spontaneous magnetization M/M s 
changes continuously with respect to the strength of J2/J1. 

Next, we calculate the local magnetization (S*^) to investigate the spin configurations in 
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FIG. 5: (Color online) Local magnetization (S?^) at each sublattice £. Panels (a) and (b) are results 
for J2/J1 =1.8 and 3.5 respectively. These results are obtained from our DMRG calculations for 
AT = 144 (z =1,2, 24). 

these two magnetic states, where (A) denotes the expectation value of the physical quantity 
A and S*g is the ^-component of S^. Figure |5] depicts our results for a system size N = 144 
on the lattice depicted in Fig. [TJb) under the open-boundary condition; Fig. EI a) and Fig. 
Mjoi) correspond to the case of J2/J1 =1.8 and 3.5 respectively; we use green inverted triangle 
for £ =A, blue pentagon for £ = A', red circle for £ = B, black cross for £ =C, aqua triangle 
for £ = C, and purple square for £ =D. In Fig. G^a), we find clearly incommensurate 
modulations with long-distance periodicity in the behavior of the local magnetization at 
the -B-sublattice sites. Therefore, we conclude that the intermediate-magnetization phase 
of < M/M s < 1/3 is the NLM ferrimagnetic one. On the other hand, in Fig. 0(b), 
we observe the uniform behavior of upward-direction spins at sublattice-sites B, C, C and 
D and downward-direction spins at sublattice A and A' as illustrated in Fig. H^b). The 
spin configuration in the phase of M/M s = 1/3 can be understood from the viewpoint of the 
MLM theorem because the present strip model corresponds to the 5 = 1/2 antiferromagnetic 
Heisenberg model on the diamond chain in the limit of Jij J\ = 00. Therefore, it is naturally 
lead to the conclusion that the phase of M/M s = 1/3 is the LM ferrimagnetic one. 

Finally, we discuss the relationship between the ground-state properties of the original 
Union Jack model depicted Fig. \T( a) and those of the strip model depicted in Fig. [TJ(b). 
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It is confirmed that the schematic spin configuration in the Neel-like spin liquid state de- 
picted in Fig. Hta) is consistent to that in the Neel state depicted in Fig. [2j We also 
confirm that the schematic spin configuration depicted in Fig. H(b) agrees completely with 
that in the semistriped-ferrimagnetic state of the original Union Jack model as shown in 
Fig. In addition, there exists the intermediate-magnetization state where the sponta- 
neous magnetization is gradually changed in the ground-state of the both models although 
the incommensurate modulation with long-distance periodicity has not been confirmed in 
the case of 2D Union Jack model. Therefore, one finds that the the ground-state phase 
diagram of the S = 1/2 antiferromagnetic Heisenberg model on the Union Jack strip lat- 
tice is qualitatively consistent to that of the S = 1/2 antiferromagnetic Heisenberg model 
on the original Union Jack lattice. The intermediate-magnetization state of the original 
Union Jack lattice have been understood from the view point of the classical configuration, 
namely " canted-state" . However, our numerical results of the Union Jack strip model leads 
to the possibility that the intermediate-state of the original Union Jack lattice is also the 
NLM ferrimagnetic one whose characteristic behavior in the local magnetization originates 
from pure quantum effects. The future studies are desirable to confirm the presence of the 
incommensurate modulation in the canted-state of the 2D Union Jack model. 



V. CONCLUSIONS 



We have studied the ground-state properties of the S = 1/2 antiferromagnetic Heisenberg 
model on the Union Jack strip lattice depicted in Fig. Hfb) by the ED and DMRG meth- 
ods. Our numerical calculations have revealed that the intermediate-magnetization state 
occurs between the Neel-like spin liquid state corresponding to the Neel state of the original 
2D model and the LM ferrimagnetic state which agrees with the semistriped-ferrimagnetic 
state of the original 2D model. In this intermediate-magnetization state of this strip model, 
the spontaneous magnetization changes gradually with respect to the strength of the inner 
interaction. We have also found the existence of the incommensurate modulation with long- 
distance periodicity of the local magnetization. From the finds of these characteristic be- 
havior, it has concluded that the intermediate state of this strip model is the NLM ferrimag- 
netism. These results naturally lead to the expectation that the intermediate-magnetization 
state of the original model is also the NLM ferrimagnetic one. 
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